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Abstract 

In the simplest little Higgs model the new flavor-changing interactions between heavy neutrinos 
and the Standard Model leptons can generate contributions to some lepton flavor violating decays 
of Z-boson at one-loop level, such as Z — )• r^/i^, Z — t- r^e^, and Z — t- /i^e^. We examine the 
decay modes, and flnd that the branching ratios can reach 10^'' for the three decays, which should 
be accessible at the GigaZ option of the ILC. 
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I. INTRODUCTION 



Little Higgs theory [l[ has been proposed as an interesting solution to the hierarchy prob- 
em. So far various realizations of the little Higgs symmetry structure have been proposed 
5|, which can be categorized generally into two classes js]. One class use the product 
group, represented by the littlest Higgs model jst, in which the SM SU{2)l gauge group is 
from the diagonal breaking of two (or more) gauge groups. The other class use the simple 
group, represented by the simplest little Higgs model (SLHM) j4|, in which a single larger 
gauge group is broken down to the SM SU{2)l. The flavor sector of little Higgs models 
based on product groups, notably the littlest Higgs model with T-parity (LHT) 5^, has 
been extensi vely studied i?!. Recently, some attentions have been paid to the flavor sector 
of SLHM Is-lO]. 



The lepton flavor violating (LFV) decays of Z-boson can be a sensitive probe for new 
physics because they are extremely suppressed in the SM but can be greatly enhanced in 



new physics models llNl3|. The experimental limits obtained at LEP 1J| are 



BR{Z T^fi^) < 1.2 X 10"^ 
BR{Z r^e^) < 9.8 x 10"^ 
BR{Z n^e^) < 1.7 X 10"^ 



(1) 



The next generation Z factory can be realized in the GigaZ option of the International 
Linear Collider (ILC) 15| . About 2 x 10^ Z events can be generated in an operational year 
of lO^s of GigaZ. Thus the expected sensitivity of GigaZ to the LFV decays of Z-boson 
could reach jl6 | 



BR{Z r^fi^) ~ K X 2.2 X 10"^ 
BR{Z r^e^) ~ k x 6.5 x 10^^ 
BR{Z -> /i^e^) ~ 2.0 X 10"^ 



(2) 



with the factor k, ranging from 0.2 to 1.0. Therefore, GigaZ can offer an important oppor- 
tunity to probe the new physics via the LFV decays of Z-boson. 

The SLHM predicts the existence of heavy neutrinos, which have flavor- changing cou- 
plings with the SM leptons mediated respectively by the SM gauge boson and the new 
heavy gauge boson X"^. These couplings can give great contributions to Z-boson decays 



Z — )■ T /i^, Z —7- r e^, and Z — )■ /i at one-loop level. In this paper, we will calculate 
the branching ratios of these decay modes, and compare the results with the sensitivity of 
GigaZ and the present experimental bounds, respectively. 

This work is organized as follows. In Sec. II we recapitulate the SLHM. In Sec. Ill we 
study respectively the decays Z — t- t^/x^, Z — t- r^e^ and Z — )■ /i^e^. Finally, we give our 
conclusion in Sec. IV. 



II. SIMPLEST LITTLE HIGGS MODEL 



0. 



The SLHM is based on \S\J{k>) x ?7(1)x] global symmetry [4J. The gauge symmetry 
^[/(S) X f/(l)x is broken down to the SM electroweak gauge group by two copies of scalar 
fields $1 and $2? which are triplets under the ^[/(S) with aligned VEVs /i and j^- The 
uneaten five pseudo-Goldstone bosons can be parameterized as 



$1 



where 



e 



1 

7 








e 







(3) 



E 



\ 




/ft 



+ 



1 
1 
1 



(4) 



/ = \/ /i + /I = tan/3 = /2//1. Under the SU{2)l SM gauge group, 77 is a real 

scalar, while H transforms as a doublet and can be identified as the SM Higgs doublet. The 
kinetic term in the non-linear sigma model is 

2 



i=i,2 



(5) 



where Qx = gtw/ a/1 — and = tan 6w with 6w being the electroweak mixing angle. 

As $1 and $2 develop their VEVs, the new heavy gauge bosons Z', Y^, Y^"^ and get 
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their masses after eating five Goldstone bosons, 

gf f, 



Mx = ^{l 



V2 \ 4/^ 



v/3^V 16/2 

My = ^. (6) 

The gauged SU (3) symmetry promotes the SM fermion doublets into SU (3) triplets. For 
each generation of lepton, a heavy neutrino is added, whose mass is 

rriNi = fs^X^- (7) 

Where i = 1, 2, 3 is the generation index and is the Yukawa coupling constant. 

After the EWSB the light and the heavy neutrino of the same family have the mixing, 
which is parameterized by 5„ = — . The mixing angel 6^ is experimentally constrained 
to be small jl7|, and taken as a typical upper limit 6^ < 0.05 following the ref. jsj. Besides, 
there is family mixing as long as the Yukawa matrix of heavy neutrinos and that of leptons 
are not aligned. This can induce the lepton fiavor-changing interactions of charged currents 
proportional to V^^'' NLi'-j^X^^iij and dyVi^'' NLi'l^W^^^Lj■, where V^*"' is the mixing matrix 



III. THE LFV DECAYS Z t^^jl^, Z r±e=F AND Z ^^e^ 

In the SLHM, the Feynman diagrams for Z — /i^e^ can be depicted by the Fig. 1, and 
the diagrams for Z — ?■ r^/i^, Z — )■ r^e^ are same as Fig.l, but replacing ^ and e with the 
corresponding final particles. For the 't Hooft-Feynman gauge, the fiavor-changing inter- 
actions between the heavy neutrino and lepton, mediated by the gauge bosons (Goldstone 
bosons) (x^) and (0^), can contribute to these decays. The relevant Feynman rules 
can be found in j^. 

The calculations of the loop diagrams in Fig. 1 are straightforward. Each loop diagram 
is composed of some scalar loop functions [isl which are calculated by using Loop Tools [l9 |. 
The analytic expressions from our calculation are presented in Appendix A. 



The SM input parameters relevant in our study are taken as ref. 20|. The free SLHM 



parameters involved are /, tp, the heavy neutrino mass mAr. {i = 1,2,3), and the mixing 




FIG. 1: Feynman diagrams for Z — > \x^e in the SLHM. 



matrix Vi which can be parameterized with standard form. To simply our calculations, we 
take the parameters 21 1 

1 



Si2 = VO, Si3 = VO.03, S23 = (^13 



65°, 



(8) 



which is consistent with the experimental constraints on the PMNS matrix 22|], and 8\z 
is taken to be equal to the CKM phase. To satisfy the present experimental bounds of 
i?r(yU — )■ 67) and i?r(/i — > eee), the mass splitting of the first and the second heavy neutrinos 
must be very small [8]. So in this paper we will take mTVi = vn^^ = mi = 400 GeV and 
m^Tg = in the range of 500 GeV-3000 GeV. Ref. j4| shows that the LEP-11 data requires 
/ > 2 TeV. In our numerical calculation we will take several values of for / = 2 TeV, 
f = 4: TeV and / = 5.6 TeV, respectively. 

In Fig. |2l Fig. [3] and Fig. HI we plot the decay branching ratios of Z ^ t^/^^, ^ — t- r^e^ 
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FIG. 2: The branching ratios of Z — )■ r^/n^ versus 
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FIG. 3: The branching ratios of Z — )• r versus ms. 

and Z — )■ jji^e^ versus for / = 2 TeV, / = 4 TeV and / = 5.6 TeV, respectively. We find 
that the branching ratios increase with the mass of the third generation heavy neutrino. The 
reason is that the decays are enhanced by the large mass splitting — mi, which increases 
as m3 gets large since we have fixed the value of mi. Besides, the branching ratios drop 
as the scale / or get large, and the reason is that the lepton flavor- changing couplings 
NLil^W^^^Lj and NiCp^ij are proportional to 6^, = -7^- 

Fig. 121 Fig. [3] and Fig. H] show the branching ratios of Z — t- r^fi^, Z — > r^e^ and 
Z — )■ /i^e^ are below the present experimental upper bounds, respectively. However, the 
ratios can be enhanced to reach the sensitivity of the GigaZ. For / = 2 TeV, = 4 and 
m3 = 2 TeV, the branching ratios can reach 10"'^ for Z — > r^fJ^^, Z — > r^e^ and Z — > /U^e^, 
which exceed much the sensitivity of GigaZ. In the LHT, all the three ratios can reach 10~^ 
12| . Therefore, the LFV decays of Z-boson may be accessible at GigaZ, and thus may serve 
as a probe of the little Higgs models. 
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IV. CONCLUSION 



In the framework of the simplest httle Higgs model, we studied the LFV decays 
Z — )■ r^/i^, Z — 7- r^e^ and Z — )■ fi^e^. In the parameter space allowed by current experi- 
ments, the branching ratios of the three decays can exceed respectively much the sensitivity 
of GigaZ, which should be accessible at the GigaZ option of the ILC. Therefore, the mea- 
surement of these rare decays at the GigaZ may serve as a probe of the simplest little Higgs 
model. 
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Appendix A: The effective coupling of Z^+e 

Here we take the effective coupling of Zfi'^e~ for example. The other two couplings 
Zr^fi^ and Zr^e" can be obtained via some corresponding replacement of the analytic 
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expressions for Z[i^e . The effective coupling of Zjji^e is given by 

^VFiF2 

+T'i.yy[N,,X{W),X{W)]+T'^SF,F2H<t>),N,,v,]+T'^SF,F2[^^^^ 
+T%sv[N^,x{^),X{W)]+T-,,,f^^,^[X{W),N,]+T-,,^^^^^^ 

+r^,,^(^) TV,] + r:,,^(„3 [x{4>), TV,] , (ai) 

where the particles in the square brackets represent the particles which contribute to the 
vertex, and ^%if(j^_n) correspond to the vertexes in Fig. l{k — n). The self-energy and vertex 
contributions in the above equation are given by 

^VF.F^ = Y^[(di4^L + ciZ{PR)(-2C,p7"7"7''-27")(c2PL + d2Pii) 

-2(^e + 4t.)l'"Cp-i^{ciC2Z[PL + did2Z{^PR) + 4mF,ic2diZ{^PL 
+Cid2Z{PR)Ca + AmFAc2dlZ[PL + Cid2Z{^PR)Ca + 2mF,Co{diZ[PL 

+ciZ^PR)(2(ge + - mF2l°'){c2PL + d2PR)\ {qe, q^i, niFi, mv , mF2) , (A2) 
= '-^fidiPL + ciPR){-4Ca07^ + r-'^Ci,j^iqe + qf.r + 2{4mF-24^)Ca 
+{4.mF - 24^){qe + q^^Co - [C^pg"" - \]r - C/37^W^ + "^f)7" - ^zCpl^l'' 

+'l°'{<iix + mF){^z -ie- 4ix)Co}{c2PL + d2PR){q^i,qe,mv,mF,mv), (A3) 
^SF^F2 = Y^[t^ap7"7V(a26i4^L + ai62^/PR) + ^7"(a26i4^L + ai62^/Pil) 
+C^7^7"(aiZ{PL + biZl^PR){ie + + mF2){a2PL + ^sPi?) 
+mF,^''{bib2Z[PL + aia2Z{^PR)Cf^jf' + mF,l''{hiZiPL 

+aiZ{^PR){4e + 4iJ. + mF2){a2PL + &2f'i?)Co] {qe, q/j., mFi,ms, mF2), (A4) 
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rfss = -'-^{-'^C^pl^ia2biPL+aib2PB)-{qe + qf.rCp-f^ia2biPL + aib2PR) 
+ [-2Ca- (qe + q^.rCo] ^^(0261^^ + ai^a^i?) + mj. [ - 2C7a 
-{qe + qfj.)°'Co]{aia2PL + bib2PR)}iqij.,qe,ms,mF,ms), (A5) 

nvs = -'-^r{ciPL + diPn)[C^j^ + {4^. + mF)Co] 

x{a2PL + b2PR){qfj,,qe,ms,mF,mv), (A6) 

nsv = '-^{a^PL + b^Pn)[C^7'' + {4^ + mF)Co]r 

x(c2Pl + d2PR){qij.,qe,mv,mF,ms), (A7) 
r^ez/w = - I6^2,^(g2 _ ^2) 7" [(-^ + s'w)Pl + .^Pi^] + me) [(2i?,7^ + (2^0 

,mv,mF), (A8) 

r?eW) = - i6vr2cvt.(p2 - ml) ^^'^^^^^ + ~ ^)^e){ciC2PL + did2Pi^) - 2m^(2ijo 

,mv,mF), (A9) 

+i^iBQ){a2biPL + ai62-PR,) + mFBo{aia2PL + bib2PR)] {q^,ms,mF), (AlO) 

Tsez/H = levr^cvKb^ - m2) [(-^/^T^ + ^ei?o)(fl2biPL + aib2PR) + mFBo{aia2PL 

+bib2PR)\ {ie + m;,)7'' [(-^ + Sh/)Pl + (Pe, ^s, m^), (All) 

where = —p^, Qe = —Pe. and Pl,r = (1 =F 75)/2. The functions B and C are 2- and 
3-point Feynman integrals [19|, and their functional dependence is indicated in the bracket 
following them. The tensor loop functions can be expanded as the scalar functions 19j. In 
our calculation the contraction of Lorentz indices is performed numerically. The parameters 
appearing above are from 

Vef : zYiciPi + diPR), Vffi : tY{c2PL + ^^2^/?), 

Sef : OiPl + BiPr, Sffi : a2PL + &2-Pfi, 

ZS+S- : -P^-), ZV+S- : gvvsg^", 

ZpV;V- : -^gyyy[{p^^-p^_Yg^'^ + (p^-p^^yg>^P + {p^^-p^)^g 

Zfj2 : triZiPL + Zj,Pn), 

where V represents gauge bosons and S represents scalar particles. These couplings represent 
the seven different classes of vertices involved in our calculation. In each class of vertices, 
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the parameters ai, bi, a2, 62, Ci, di, C2, d2, gvss, gvvs, gvvv, Z[ and take different 
values for^ifferent concrete coupling. The analytic expressions of these parameters can be 
found in 
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